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We consider two-dimensional dispersions of droplets of isotropic phase in a liquid with an XY-like 
order parameter, tilt, nematic, and hexatic symmetries being included. Strong anchoring boundary 
conditions are assumed. Textures for a single droplet and a pair of droplets are calculated and 
a universal droplet-droplet pair potential is obtained. The interaction of dispersed droplets via 
the ordered phase is attractive at large distances and repulsive at short distances, which results in 
a well defined preferred separation for two droplets and topological stabilization of the emulsion. 
This interaction also drives self-assembly into chains. Preferred separations and energy barriers to 
coalescence are calculated, and the effects of thermal fluctuations and film thickness are discussed. 
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I. INTRODUCTION 

Emulsions are an important noncquilibrium state of 
bulk matter, widely used in the industry and of intrinsic 
physical interest. Although two-dimensional emulsions 
are less common, they may mediate certain membrane 
bounded process in living cells [l|, H, d, H|, control the 
kinetics of phase transitions in thin films, and shed light 
on the corresponding phenomena in three dimensions. 

One of the main questions in studies of emulsions from 
both a theoretical and experimental point of view is sta- 
bilization. Since emulsions are out of equilibrium, they 
tend to phase separate. Therefore for an emulsion to be 
useful it has to be stabilized by kinetic barriers to coales- 
cence. Typical strategies include surfactant or grafted 
polymer stabilization and could involve Coulombic or 
steric repulsion of dispersed droplets [1, @. Recently 
there has been substantial interest in topological stabi- 
lization, which may arise when droplets of disordered liq- 
uid are introduced in an anisotro pic host fluid |7 l . |q . l9|.|l0l. 

OS El, El, El EE El, EE El, El, 11, M m, HI El, HI. 

Anisotropic liquids exist both in two and three dimen- 
sions. Tilt and bond orientational order are quite com- 
mon in liquid monolayers (26j , and there is some evidence 
that hydrated lyotropic bilayers of DMPC (dimyristol- 
phophatidylcholine, saturated fatty acids chains) have 
bond orientational order as well [27], Hl[ . Quite recently, 
a new hexatic phase, with six-fold orientational order, 
has been discovered in freestanding lyotropic thin films 
as few as four bilayers thick [2^| . 

Langmuir-Blodgctt monolayers at air-water interfaces 
have long lived domains due to long-range electric dipole 
interactions between lipids at the asymmetric air-water 
interface [II]. However, these dipole interactions arc 
much reduced for bilayer systems in water, and are un- 
likely to account for the arrested phase separation ob- 
served in, e.g., Refs. [1, Hl[. Komura et. al. have con- 



*Electronic address: papcrs.korolev@gmail.com 



structcd a mean field theory which illuminates a num- 
ber of aspects of lateral phase separation in mixtures of 
lipid and cholesterol in bilayers [32| ■ The "lipid ordered" 
and gel phases are described by a simple Ising-type or- 
der parameter without, however, explicit regard for the 
orientational order that may be responsible for the slow 
coalescence. 

Two-dimensional topological emulsions have also been 
observed experimentally in thin films of smectic- C/C* 
liquid crystal [Tl|, El, El, EH- In such systems the or- 
der parameter is the projection of the tilted director on 
the plane of the film and inclusions are typically disklikc 
"bubbles" of less ordered phase (nematic or cholcsteric) 
or domains with a different number of smectic layers. 
Under appropriate conditions, these inclusions are ac- 
companied either by a topological defect in the bulk of 
the host fluid [H|,G3 or by several surface defects [H,E3|- 
The distortions in the texture lead to long-range attrac- 
tion and short-range repulsion between bubbles, and they 
drive self-assembly into chainlike and latticelike struc- 
tures [ll], El, El, E3] • Thesephenomena have been stud- 
ied theoretically before [1, El, EE EE El- Analytical 
solutions for the texture around a single droplet have 
been calculated in Refs. QI, El, HH and for the far-field 
interactions in Refs. [1, Ell- Textures around several in- 
clusions with surface defects have been obtained numer- 
ically by Bohlcy and Stannarius [n|. For bulk defects, 
short range repulsion and stability have been discussed 
in Refs. [ll, Ell- 
in this paper we discuss how orientational order and 
point defects give rise to the rich behavior seen in the 
experiments. We first review the basic properties of 
topological emulsions in Sec. [Til In the next three sec- 
tions we verify, extend and correct the results from 
Refs. d, El, E3|- In particular, wc consider droplets 
of different sizes embedded in a two-dimensional liquid 
crystal with an arbitrary orientational order parameter. 
In Sec. IIII1 we calculate textures induced by bubbles 
of disordered phase in a liquid with a generalized XY- 
like order parameter (tilt, nematic, and hexatic order 
are special cases; the relatively unexplored case of 4- 
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fold "tetratic" order is included as well.). In Sec. IIV[ wc 
obtain an analytical expression for far-field interactions 
between bubbles mediated by the aforementioned order 
parameters. We then discuss numerical computation of 
an effective droplet-droplet pair potential in Sec. [V] and 
show that this interaction leads to bubble stabilization 
in the presence of the "topological quasi-long range or- 
der" [33[ characteristic of two-dimensional (2D) systems 
with an XY order parameter. In Sec. IVII we use our nu- 
merical approach to address chaining. Finally in Sec. lVIII 
we discuss how our results are affected by thermal fluc- 
tuations and film thickness. 



II. TOPOLOGICAL EMULSIONS 



Phenomena of interest in this paper can arise in a va- 
riety of physical systems @, i, M, El 0, GJ, M, HI [H| . 
In this section we outline a model that captures many 
essential features and has a virtue of simplicity. 

A natural way to create a two-dimensional emulsion is 
to rapidly quench the system below a mixing-demixing 
transition. The defect-mediated cmulsification of inter- 
est to us here arises when the continuous phase has or- 
der parameter with XY symmetry. The system with 
both mixing-demixing and XY ordering transitions can 
be described with the following coarse-grained lattice gas 
Hamiltonian 1341. 



I j j ^ 



k B T 2 



(1) 



where K > and J > describe interactions, /i is a 
chemical potential, i and j index lattice sites, means 
a sum over nearest neighbors, ti — denotes a patch of 
disordered phase and ti = 1 denotes an ordered phase 
region. The continuous AY-like symmetry is embodied 
in the complex p-fold symmetric order parameter ipj = 
e ip8j . 

At low K (high temperature) entropy dominates and 
the system is mixed, while at high K (low tempera- 
ture) it phase separates. In addition to this Ising-type 
term the Hamiltonian includes Ay-like interactions be- 
tween sites that are occupied by the ordered phase. 
Here p = 1,2, and 6 corresponds to tilt, nematic, and 
hexatic order parameters respectively; tetratic order is 
described by p — 4. Fig. [1] shows a phase diagram ob- 
tained via real space renormalization methods for this 
Hamiltonian when K = J 34]. Thermal fluctuations 
drive the ordering transition for the continuous symme- 
try degrees of freedom {ipj} well below the Ising-type 
critical point, and also lead to the broad flat top seen 
in experiments 0, H[ and typical of the small Ising-type 
coexistence curve exponent f3 = 1/8. Because thermal 
fluctuations are strong in two dimensions, the thin sliver 
of p-fold XY ordered regime on the left side only ex- 
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FIG. 1: Phase diagram for a mixture of ordered and disor- 
dered phases. 



hibits "quasi-long-range order," with algebraic decay of 
correlations in ipi [33|. 

Well below both phase transitions the free energy of the 
ordered phase can be described by a gradient expansion 
in the angle field 6j . For p = 1 , 2 we make the one elastic 
constant approximation, 



F 



d 2 x{vef 



(2) 



This representation of the low energy physics is exact 
for p > 2. Here g is the stiffness constant, and the an- 
gle 8 is the orientational order parameter. Functional 
minimization of Eq. ([2]) leads to a simple Eulcr-Lagrange 
equation for 9, 



V 2 9 = 0. 



(3) 



One important requirement for topological stabiliza- 
tion of emulsions is the anchoring boundary energy, 
which depends on the angle between the director and 
local normal to the boundary of a bubble. Here we as- 
sume the limit of strong boundary conditions, when this 
an gle i s chosen to minimize the boundary energy 0, d, S 
Hoi. Il3l| . Tangential and homotropic boundary conditions 
are two examples when the angle is 7r/2 and respec- 
tively. Strong anchoring implies that the order param- 
eter uniformly changes by 2ir along a bubble perimeter. 
For simplicity we adopt homotropic boundary conditions 
but the results are valid for all other anchoring angles, 
because the textures only differ by a rotation of the order 
parameter by an anchoring angle at every point in space; 
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this rotation does not affect the energy of the system in 
the one elastic constant approximation of Eq. @. 

Boundary conditions at the outer boundary of the 
emulsion are also very important. We assume free bound- 
ary conditions, which simplify the calculations and are 
appropriate when there is no anchoring at the outer 
boundary. In the experiments with smectic films, the 
director at infinity is not free but is usually forced to 
point in some particular direction [ll|, E 0, 03 ■ Al- 
though it would be straightforward to repeat our calcu- 
lations with these boundary conditions, the difference is 
often unimportant because, when topological defects are 
included, both boundary conditions lead to uniform di- 
rector field at infinity and, therefore, to the same solution 
up to a rotation relative to the outer boundary. However, 
uniform boundary conditions at infinity must be used if 
one studies the dependence of interaction energy on the 
angle between the line connecting the droplets and the 
orientation of the director at the outer boundary. The 
free boundary conditions correspond to the angle that 
minimizes the energy of the system. 



fects we have only to determine the distance R between 
the center of the droplet and the defects. On dimensional 
grounds, R = Sa, where S is dimensionless. To calculate S 
we have to find a solution of Eq. that satisfies the 
boundary condition and has appropriate singularities for 
the pattern of defects specified above. This is most easily 
accomplished by the method of images. The solution is 



6 = 2(f) Q - - V(> ? +<#•) + const, (4) 

y i=l 

where the {<pi} are polar angles around the defects and 
the {(/>[} are polar angles around the images obtained by 
inversion in the circle. Technical details of the derivation 
are discussed in Appendix [SJ The inverse problem of 
textures inside a circular droplet of ordered phase em- 
bedded in a disordered isotropic phase is discussed in 
Appendix [B] 

From Eqs. (0| and ^ we obtain the free energy as a 
function of R, namely, 



III. SINGLE DROPLET 

First we consider a single circular droplet of isotropic 
phase with radius a embedded in a two-dimensional liq- 
uid crystal with p-fold symmetric XY order. We assume 
a strong line tension preserves the circular shape. In 
the absence of topological defects, the only solution sat- 
isfying Eq. ^ and the boundary condition is 9 = </> , 
where 4>q is the polar angle relative to the center of 
the droplet. Thus, from the outside, the droplet looks 
like a point defect with charge +1. From Eq. ([2]) wc 
find that the free energy of this texture is proportional 
to irgh\(L / a), where L is the system size. However, the 
free energy can be lowered significantly below this loga- 
rithmic divergence if one allows for an appropriate num- 
ber of point defects, so that the net charge of the droplet 
and the defects is zero. 

The order parameter we are considering can be rep- 
resented by an arrow with directions differing by 27r/p 
identified. This order parameter allows for point defects 
with topological charge n/p, where n is an integer, which 
means that the order parameter rotates by 2irn/p on any 
contour surrounding the defect. To make the free energy 
of the system finite in the presence of a droplet, we need 
the charges of the defects to sum up to —1. Upon noting 
that the core energy scales as n 2 and that many defects 
with smaller charges screen better than few defects with 
larger charges, we conclude that n = 1 and study p de- 
fects with charge —1/p. From the symmetry considera- 
tions and because defects repel each other it is clear that 
they must form a regular p-gon around the droplet, at 
least for p of physical interest, p < 6. This behavior was 
also found numerically by Fukuda and Yokoyama 
for p — 2 considered there. 

To completely specify equilibrium locations of the de- 
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where c is the core radius and E c is the core energy. Upon 
minimization, Eq. (J5]) gives (with S = R/a) 



3p + l 
p+l 



(6) 



and a free energy that no longer diverges with the system 
size L, 



F=^hJ 



S 3pA 



p \p(S 2 P - 1) c 



pE c 



(7) 



Our results agree with the calculation in Ref. [9|] and 
experiments by Cluzeau et al. which apply to the 
case p = 1 and with the calculation in Ref. [lf|, which 
apply to the case p = 2; they are also consistent with 
the numerical work discussed in Sec. |V] In another ex- 
perimental system with p = 1, a different texture has 
been observed [H, [3| ■ Instead of one topological defect 
in the bulk, a droplet has two and sometimes four sur- 
face defects of strength 1/2 and 1/4 respectively. In the 
case of the one elastic constant approximation and free 
boundary conditions at infinity considered here, such de- 
fects are located equidistantly along the circle because 
they repel each other. The texture could be obtained by 
the method discussed in Appendix [Al in this case the de- 
fects and images are located at the same points. For two 
half-defects the elastic energy is decreased by 31n(2)7rg 
compared to the texture with one defect outside the bub- 
ble, but the anchoring energy is increased because the 
boundary conditions must be violated around two singu- 
lar points. The energy difference between the alternative 
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textures does not depend on the size of the droplet and 
the winning configuration is determined by the competi- 
tion between elastic and anchoring energies. We assume 
here anchoring energies large enough to exclude surface 
defects on the bubble. 

In summary, we have obtained the exact texture 
around a single droplet with strong anchoring boundary 
conditions on the surface of the bubble and free boundary 
conditions at infinity. One important conclusion is that, 
in addition to line tension and bulk contributions to the 
droplet free energy, there is a new size-dependent con- 
tribution to the droplet free energy (irg/p)Yi\(a/c). This 
additional energy will alter the usual estimates for nuclc- 
ation and growth of droplets of isotropic phase in terms 
of bulk and surface energies. 

It is also important to emphasize that the penalty for 
violating the strong anchoring boundary condition on the 
droplet itself, by superimposing, say, a uniform texture 
around the bubble, grows linearly with a. However, the 
creation of several topological defects only bears the fixed 
cost of the core energies. Thus, small droplets may not 
be accompanied by the defects, while large ones are. 

While we have shown that a single droplet acquires 
companion defects, demonstrating topological stabiliza- 
tion of emulsions requires more. We now consider the in- 
teraction of two such domains as a function of the center- 
to-center distance r between them, and show that it leads 
to large barriers to coalescence. 



IV. FAR-FIELD INTERACTIONS 
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FIG. 2: Geometry of the problem in the far field limit, r S> 
01, a2 • Droplets of radius a\ and a 2 are surrounded by regular 
p-gons of defects, and the pair potential depends only on their 
orientations and the separation r. 

ued functions or to introduce many branch cuts, both of 
which complicate the problem. An alternative approach 
is to introduce auxiliary single valued fields tp and \ such 
that, 

V<0 = ViV - ea Vj-x, (8) 

where is the antisymmetric tensor in two dimensions. 
It is possible to decouple the auxiliary fields by demand- 
ing that ip is constant on the surface of the droplets 
(the constants do not have to be the same on different 
droplets), which implies that V?/> is normal to the bound- 
aries. The energy is then given by 



Similar to our analysis of a single droplet, we expect 
to have 2p defects with charges —1/p to neutralize the 
topological charge of two droplets with strong anchoring 
boundary conditions. Unlike the single droplet case, we 
cannot predict the location of the defects based only on 
the symmetry; we need other methods to find equilib- 
rium positions of the singularities for each droplet sepa- 
ration r to compute an effective droplet-droplet pair po- 
tential V(r). 

In the far- field limit, when the distance between the 
droplets is much greater than their sizes, one can still as- 
sume that the defects form a regular p-gon around each 
of the droplets; the sizes of these p-gons remain unaf- 
fected, to leading order in a/r. What one cannot neglect 
is the relative orientation of the p-gons. Due to the ro- 
tational symmetry of an isolated droplet, its energy does 
not depend on relative orientation of the p-gon and the 
droplet. For two droplets the symmetry is reduced, and 
the interaction should set the preferred orientations of 
the p-gons. The geometry of the problem is summarized 
in Fig. H 

Thus, we have to solve Eq. ((3]) subject to the homo- 
topic boundary conditions used for a single droplet and 
specified locations of the point defects. The calculation is 
not straightforward because 9 is defined only up to 2tt/p 
rotations. This leads one either to consider multival- 



F= 9 -JmfSx+ 9 - J(V X ) 2 d 2 x, (9) 
where the cross term vanishes since 

- / (€«,■ V^V jX )d 2 x = / { X e tJ V l V^)d 2 x+ 

J \ (10) 

/ (x^jUjVi^dl. 

The first term on the right hand side is zero due to 
antisymmetry of e^ . The second term vanishes for a 
similar reason, because we choose 'Vtp to be parallel to 
the outward unit normal to the boundary n. We also 
assume that tp and \ decay fast enough that there is no 
contribution at infinity. 

It follows from Eq. ([9]) that both ijj and \ satisfy 
Laplace's equation. One can also show that topologi- 
cal defects in the 0— field correspond to quantized sin- 
gularities in \~ field and that for any texture there al- 
ways exist ip and x which satisfy Eq. ([8]). The approach 
of dividing the gradient of the order parameter in two 
fields is similar to the electromagnetic analogy exploited 
in Ref. @. 

To find the interaction of droplets we solve the follow- 
ing set of equations: 
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V 2 V> = 0, 



V 2 x=-£^(x- di ), 

n — * 



2p 



P 



(11) 



(12) 



where x labels points in space and di refers to the posi- 
tions of the defects. We choose boundary condition for \ 
on the surface of the first droplet to be 
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2.190 
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2.167 


0.996 
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154 


2.11 


0.77 



TABLE I: Values of A{p) for Eq. (Tl6]) . preferred separa- 
tions r m i n and barriers to coalescence A for physically rel- 
evant values of p and a\ — a^. Values of r m i n /a and A/g are 
accurate up to the last digit. 



Vx-n = l/a x . 



(13) 



This condition guarantees that the order parameter uni- 
formly changes along the boundary by 2tt. However, the 
director generally forms a nonzero angle ot\ with the nor- 
mal. This discrepancy can be easily fixed by imposing a 
boundary condition on Eq. (jll[) for tp that effectively ro- 
tates the order parameter to match the normal as follows: 



= -ai(r, 0i, 2 , ai, a 2 ). 



(14) 



We impose analogous boundary conditions on the surface 
of the second droplet. 

The equations for tp are solved exactly in Appendix O 
in the far-field limit, the tf> contribution to the interac- 
tion energy [i.e. the contribution from the first term 
of Eq. ([9])] reduces to 



irg(a 2 - aif 



In 



(15) 



where a\ and a 2 are functions of r, a\ , a 2 and positions of 
the defects, which are specified by 0i and 02 in the far- 
field limit. Equation (|15p has strong l/ln(r) dependence, 
but as we will see later one can always arrange the defects 
in a way to set u\ = a 2 and eliminate V^. 

For the \ field; Eqs. (JT5J) and (JTSJ) could be solved 
approximately by the method of images. The number 
of images is infinite, but only a few of them contribute 
to the leading order. Apart from the p images inside 
each droplet that we considered in the previous section, 
we have to include images inside the second droplet of 
the whole structure in and around the first droplet (de- 
fects, their images, and +2 charge at the center of the 
first droplet; see Appendix |A")) . as well as images inside 
the first droplet of the whole structure near the second 
droplet. 1 We use this solution to compute the energy, 
which we then Taylor expand to the leading order in 



1 Unlike Ref. @, we find that one has to consider more than 2p 
images 



powers of 1 /R. The result for the x P ar t of the droplet- 
droplet pair potential is 



V x (r^ 1 ,^ 2 ) x =(-iy7rA(p)g 



5p 



x cos[p(0i + 2 )] 



r 2p 



(16) 



where A{p) is given in Table [IJ 

There are two important features of these far-field pair 
potentials. First, for any value of p there exist an- 
gles 0i and 02 such that V x (r, 0i, 02) < 0, and there- 
fore V x (r) can lead to an attraction, unlike V^,, which 
is always repulsive. Second, there is a one-dimensional 
family of degenerate energies, because one can always 
add an arbitrary angle to 0i and subtract it from 02 
without changing V x . Similarly, depends only on the 
difference between u\ and a 2 - 

To complete the analysis we only have to 
find a± and a 2 in terms of 0i,02, and r. In prin- 
ciple this dependence is given by the set of images 
described above but the explicit expressions are some- 
what complicated. It is, however, easy to see that 
when ai = a 2 and 0i = 02 mod 27r/p the contribution 
to the energy from the ip field should be small due 
to symmetry Moreover, one can show analytically 
that for this orientation of the p-gons, falls of at 
least as r~( 2p+1 ) / ln(r), which is much faster than V x . 
Thus, for a\ = a 2 we can minimize the total interaction 
energy V = + V x by requiring that 



mod 



[l + (-l) P ]7T 
2p 



mod 



(17) 



(18) 



where Eq. (fT7|) follows from Eq. p^]) . If p = 1 the solution 
is 0i = 02 = or 0i = 02 = 7r and if p is even 0i = 
02 = n/(2p) + 7rfc/p, where k is an integer. For different 
size droplets Eq. (JTHJ) still holds but Eq. flT7|) must be 
generalized to an implicit equation 



6 



ai(r, ai,a 2 , <t>\ t <fa) = a^r, ai, as2, ^i, <fe)- 



(19) 



Generically, Eqs. (|18|) and (fl9|) have a solution, which 
minimizes both and V x . The large r pair potential 
for this favorable relative orientation is then attractive 
and is given by 



V{r) 



-7rA(p)g \P 



1 

6p 



p p 

a l a 2 
r 2p 



(20) 



It is worth noting that V x obtained here agrees with 
Ref. Q for p — 1 considered therein, but we have a dif- 
ferent expression for V^. Wc think that our result is 
more physical because is repulsive, does not depend 
on the system size, and decays for large r, unlike the 
potential calculated in Ref. Q. Our result also gener- 
alizes the work of Tasinkevych, Silvestre, Patricio, and 
Telo da Gama [16| , who considered the case of p = 2, a\ = 
a 2 , and 4>i = 02- Since in this case = 0, the authors 
did not discuss the role of V^, in the pairwise interac- 
tion. The superposition principle does not hold for rea- 
sons related to footnote HJ as a result, the pair potential 
in Ref. flU ] differs from ours in its dependence on S. 



V. NUMERICAL CALCULATION OF 
PAIRWISE INTERACTION 

As the droplets get closer to each other and subleading 
terms become important, we expect the interaction be- 
tween the clouds of defects to give rise to a repulsive force 
somewhat similar to electron mediated repulsion between 
atoms of ordinary matter. To check this conjecture we 
must extend our calculation to r w ai,a%. To that end 
we devised an algorithm of solving Eq. (1121) subject to the 
boundary conditions (|13p and specified positions of the 
defects; we used Eq. (|C5[) from Appendix [C] to account 
for the energy of the ip field. This procedure was carried 
out both by solving Laplace's equation on a lattice and by 
iterative evaluation of a large number of images. The re- 
sults of these two methods were mutually consistent and 
the relative error of energy could be made as low as 10~ 5 . 
An output of one such calculation for p = 1 is shown in 
Fig. [31 We then used a downhill simplex method [35[ to 
minimize the free energy with respect to the positions of 
the topological defects; different starting simplexes were 
used to check that the algorithm converges to the correct 
global minimum. The procedure was repeated for several 
separations r until a reliable graph of the pair potential 
was obtained. This calculation was done for p = 1, 2, 4, 6 
and different ratios of 02/01. 

A result of one such calculation for p = 1, a 2 /ai = 1 
is shown in Fig. [4] where the inset shows locations of the 
defects at the minimum of the pair potential. At large 
separations, droplets attract and V(r) = — 9irg(a/r) 2 in 
agreement with Eq. (f2"0")) . Note, there is no logarithmic 




FIG. 3: Equilibrium texture for a 2D, p = 1, liquid crystal 
for equally sized droplets at the preferred separation, with 
free boundary conditions at infinity. A defect with topolog- 
ical charge —1 appears just to the left of each droplet. The 
brightness represents the square of sine of the order parame- 
ter, which should be similar to what one would observe with 
crossed polarizers. 



contribution from because a\ = a-i for the equilibrium 
configuration of defects. At small separations, V(r) is 
repulsive and when the droplets touch it tends to infinity, 
thus inhibiting coalescence. A distinct feature of this pair 
potential is a minimum which implies that at sufficiently 
low temperatures two droplets come together to form a 
stable dimer, as shown in Fig. [3j 

As wc showed in the previous section, the defects 
for p = 1 prefer to lie on the line connecting the centers of 
the circles when the droplets are far apart (see the inset in 
Fig. ID) From our simulation we find that this behavior 
persists for all separations, r and relative droplet sizes 
studied 02/01 (wc explored 0.1 < 02/01 < 1), because 
for this configuration of the defects ot\ — cti and is 
zero. Therefore, even when the droplets come very close 
there is a defect sitting between them, which naturally 
leads to a diverging pair potential as r — ► (ai + 02). 
This divergence will be cut off either by the core radius, 
or when our strong anchoring boundary conditions are 
relaxed. Textures and r m j n obtained here agree reason- 



13] 
17] 



ably well with measurements by Cluzeau et al. | 
for ap = 1 system and with simulations in Refs. 
for p = 1 and p = 2 systems. 

For p = 2,4, and 6 the pair potential also has a mini- 
mum, but V(r) remains finite when droplets touch. The 
difference arises because there is a centrosymmctric con- 
figuration of the defects that minimizes V x and elimi- 
nates 1/0. In addition, this configuration (see the insets 
of Fig. does not have a defect between the droplets, the 
cause of the divergence for p = 1 . Examples of such pair 
potentials are shown in Fig.[5]for p = 2, 6, and 02/01 = 1. 

We quantify barriers to coalescence by introducing a 
free energy difference A between states when droplets 
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FIG. 4: Pair potential for p = 1 and ai = a2. The dots on 
the curve are the actual data points while the line is just a 
guide to the eye. The inset shows equilibrium locations of the 
two defects at the preferred separation, i.e. at the minimum 
of the pair potential. The defects prefer to lie on the center 
to center line, as shown on the inset, for all values of r. The 
behavior for other values of 02/01 is qualitatively similar. 



touch (r = ai + 0,2) and when they are at the preferred 
separation r m - m , which provides a lower bound on the ac- 
tivation barrier, and leads to e _A / feBT -fold decrease in 
the rate of coalescence. It is important to know how A 
depends on the ratio aija\ of the two droplet sizes to un- 
derstand the dynamics of coalescence. This dependence 
is summarized in Fig. for p = 2; we find similar be- 
havior for p = 2,4, and 6. One can see that the more 
asymmetric the size ratio the higher the barriers to coa- 
lescence, which implies almost equally sized droplets are 
more likely to coalesce than droplets of different sizes. 

For all physically relevant values of p and ra- 
tios of 02/01 studied, we found barriers to coalescence 
and a minimum of the pair potential, which implies topo- 
logical stabilization and dimcrization of isolated droplet 
pairs. Tablc[T] (see Scc. HV[) summarizes preferred separa- 
tions and minimum barriers to coalescence for different 
values of p and a\ = d2- 

The formation of droplet dimcrs has a number of com- 
mon features with the chemistry of diatomic molecules. 
For equal droplet sizes the defects are divided evenly be- 
tween the droplets similarly to a nonpolar covalcnt bond. 
If the sizes are slightly different, the defects shift toward 
the smaller droplet as in a polar covalent bond. Eventu- 
ally, when the radii differ by about a factor of 2 or more, 
the smallest droplet annexes one of the defects and the 
bond becomes ionic (see FigJT]). 

Our calculations not only explain the stability of emul- 
sions, but also suggest an experimental signature of topo- 
logical emulsions. For example, hexatic order can be diffi- 
cult to detect directly, because it does not couple directly 
to the polarization of light; detection via x-ray scattering 
experiments on two-dimensional films can be challenging. 




a Center to center distance, r/a 

0.8r- 




u ~ 2 2 22 2.4 2.6 

Center to center distance, r/a 

FIG. 5: Pair potentials for (a) p — 2, ai = 02 and (b) p = 
6, ai =0,2- The dots on the curves are the actual data points 
while the lines are just guides to the eye. The insets show 
equilibrium locations of the defects at the preferred separa- 
tion, i.e. at the minimum of the pair potential. The behavior 
for other values of 02/01 is qualitatively similar. Note the 
centrosymmetric configuration of the defects. 



In contrast, dimerization is observable by ordinary light 
microscopy and our theory predicts a universal equilib- 
rium separation that depends only on the sizes of the 
droplets but not on the properties of the material. For- 
mation of dimers with separation that agrees with the 
theory presented here (there are no adjustable parame- 
ters) would be a strong indicator of the presence of p-fold 
order surrounding coexisting droplets. In addition, our 
calculation of the pair potential could be used to extract 
the stiffness constant g from experimental measurements 
of thermal fluctuations about the preferred separation. 
This signature could also be important for the order pa- 
rameters (p = 1 and p — 2) that couple directly to light, 
because precise measurements of the texture and location 
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1.3 



0.25 0.5 75 1 
a 2 /3 1 

FIG. 6: Barriers to coalescence for p — 2 as a function of the 
ratio of droplet sizes, A(a2/ai). The dots on the curve are 
the actual data points while the line is just a guide to the eye. 
Note that A(a2/oi) = A(ai/a2). 




FIG. 7: The smaller droplet annexes a defect as the separa- 
tion between the circles decreases. Both (a) and (b) show 
equilibrium configurations of the defects for two droplets 
with ai = 2d2- In (a) r = 3.3ai, and each droplet has a 
companion defect, while in (b) r = 3.2ai, and both defects 
are located near the smaller droplet. The figure shows an- 
nexing for p — 1, but similar behavior has been observed 
for p = 2, 4 and 6. 



of defects could be more difficult than measurements of 
the separation between the droplets (l3| . 



FIG. 8: Minimal energy configuration of three equally sized 
droplets, p — 2. The inter droplet spacing is 2.19a, which is 
very close to the size of a dimer. 



VI. CHAINING 



Even though understanding the interaction of two 
droplets is important, real systems can be more subtle, 
because simple pairwise interactions sometimes lead to 
complex behavior of matter as a whole. Indeed, many 
experiments report self-assembly of droplets into chains 
at low concentrations and lattices at high concentra- 
tions [ll], O O 03 ■ Chaining could be explained by 
texture- induced dipole-dipole interactions in the presence 
of an in-plane aligning field due to boundary conditions, 
similar to ferrofluids [5[. However, this explanation seems 
insufficient in two ways. First, at small separations, as 
one might expect in chains, the interaction is no longer 
dipole-dipole as we showed above. Second, the expla- 
nation appeals to an aligning field or aligning boundary 
conditions at infinity, dismissing a possibility of chaining 
in their absence. 

Unfortunately, our numerical method of images be- 
comes prohibitively complex for more than a few droplets 
and requires free boundary conditions. In addition, 
for several droplets one has to solve for ip numerically. 
Nevertheless, there are two important cases where the 
progress can be made. 

The first case is chaining for p = 1, where we know 
from experiments by Cluzeau et al. [Tl| that defects lie 
along the chain between the droplets. For this config- 
uration is zero (which probably drives the chaining) 
and we can find distance between neighbors in a chain of 
equal droplets, which is 2.89a. 

The second case is chaining for p = 2, where the tex- 
ture around three and four equal droplets that mini- 
mizes V x happens to set to zero as well. This con- 
figuration of defects is shown in Fig. [8] and is analogous 
to a polymer such as polyethylene in chemistry. In fact, 
one can view the sealing defects on the ends as free rad- 
icals that facilitate the reaction of chain polymerization. 
An important conclusion is that in this case chaining is 
driven solely by the interaction of topological defects - 
no external field is required. 
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FIG. 9: Average separation between the droplet and the de- 
fect as a function of temperature for p = 1. The separation 
approaches a value given by Eq. as T — > 0. 



VII. EFFECTS OF THERMAL FLUCTUATIONS 
AND BOUNDARY CONDITIONS 

It is well known that thermal fluctuations are espe- 
cially important in two dimensions. In this section we 
discuss how they modify the results presented so far. In 
the one elastic constant approximation, considered in this 
paper, the effects of temperature depend on only one di- 
mcnsionlcss parameter, h = ksT/g. Even though the 
temperature can be varied only in a limited range be- 
tween neighboring phase transitions, in most liquid crys- 
tals, g can be changed by about two orders of magnitude 
by varying film thickness in the experiments of Ref. pd| . 
In fact, g « UbTcN , where T c is the temperature of the 
ordering phase transition in a single layer, and iV is the 
number of layers in the film [9( . Thus one can roughly 
interpret x^ 1 as N. 

Even though thermal fluctuations lead to many impor- 
tant effects in liquids with a continuous XY symmetry, 
below the Kosterlitz-Thouless transition most of them 
can be absorbed into a renormalized stiffness g that now 
depends not only on temperature but also on the length 
scale [33[ . We neglect this dependence because it is rather 
weak, but we allow the defect positions to fluctuate in 
space according to the Boltzmann distribution. Then we 
compute thermal averages of interest either by direct nu- 
merical integration weighted by Boltzmann factors, or by 
Monte Carlo simulations. 

First we analyze a single droplet for p = 1. Due to 
the highly asymmetric shape of the potential in Fig. Q] 
we expect that on average the defect is located further 
from the droplet than Eq. (J6]) predicts. This effect is 
significant even for small values of x as shown in Fig. [5J 
and could be observable experimentally. 

Effects of thermal fluctuations are even more dramatic 




R/a 

FIG. 10: (Color online) Probability density for the separa- 
tion between a defect and the droplet for p = 6. The inset 
shows that the six defects execute random walks; different col- 
ors represent different defects. The yellowish annulus around 
the boundary arises because random walks overlap near the 
droplet, where the density of the defects is very high. The 
simulation was done at xT x = 25 with 10 6 Monte Carlo steps, 
but only each tenth step is shown for clarity. 

for higher values of p. Figure [TO] shows the probabil- 
ity density for the separation between a defect and the 
droplet when p = 6; the inset shows random walk of the 
six defects. Similar to the blurred electron wave function 
of atomic physics, topological defects at high tempera- 
tures form a "defect cloud" around the droplet rather 
than a simple hexagon. 

For two droplets, the quantity of interest is an effec- 
tive pair potential for nonzero temperatures. By effective 
potential we mean thermal average of the bare pair po- 
tential with respect to the positions of the defects. Fig- 
ure [11] shows this quantity for p = 1 and >r _1 = 4. Upon 
recalling that k~ x w N, one can sec that even for a film 
of only four layers thick there is still a minimum and 
barriers to coalescence, which suggests that our results 
are robust to thermal fluctuations. It is also interesting 
that r m i n increases with temperature (and decreases with 
film thickness). The trend is the opposite of what was 
found experimentally in Ref. (l3| . suggesting that some 
other effect must account for the observations. Similar 
results were obtained for p = 2. 



VIII. CONCLUSIONS 

We have shown that two-dimensional dispersions of 
isotropic phase in a liquid with a continuous symmetry 
XF-like order parameter are topologically stable, pro- 
vided there are strong anchoring boundary conditions 
on the droplets. We have calculated the texture for 
one and two droplets and determined the pair potential 
for p = 1,2,4, and 6. A cloud of defects surrounding each 



10 




3 4 5 

Dropet to droplet distance, r/a 



Research Science and Engineering Center through Grant 
No. DMR-0213805. 



APPENDIX A: ENERGY OF A SINGLE 
DROPLET 

In this appendix we derive Eqs. ([5]) and ([6]) for the en- 
ergy and position of the satellite defects around a droplet, 
starting from just before Eq. (T4J. To display the solution, 



we first introduce some notation. Let At 



1, 



,P be 



the positions of the defects created in the liquid crystal, 
and A\ , i = 1 , . . . , p their images in the circle C of radius a 
centered at the origin O, which represents the surface of 
the droplet. The images are obtained by inversion, i.e. 
the image (x' , y') of (x, y) in C is given by 



FIG. 11: Thermally renormalized pair potential for p = 
1, til = 02, and x" 1 — 4. The dots on the curve are the actual 
data points while the line is just a guide to the eye. Energy of 
the droplets is measured in units of g. Unlike Figs. f4] and 
we did not subtract the energy of two isolated droplets. 



droplet insures topological charge neutrality and con- 
tributes an additional term to the free energy, which de- 
pends logarithmically on the size of the droplet. At large 
separations, droplets of isotropic phase attract while at 
small ones they repel. The defect-mediated pair potential 
clearly indicates kinetic barriers to coalescence and favors 
dimcrization. We developed a fast numerical method for 
obtaining the pair potential for droplets of similar and 
disparate sizes and used it to calculate activation barri- 
ers and the sizes of dimers. The latter results could be 
used to test experimentally for the presence of "hidden" 
order parameters with, e.g., p = 4 and p = 6. We also 
analyzed chaining and showed that it could be driven 
not only by an aligning field due to the boundary condi- 
tions @, [3 but instead solely by the interaction between 
the droplets and topological defects for p = 2. Finally, 
we addressed the effects of thermal fluctuations on the 
texture around a single droplet and on the pair poten- 
tial and indicated how the preferred separation could be 
affected by the thickness of the film. 
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At each position Ai we define an angular variable <f>i(X), 

> 

which is the angle between axis OAi and vector AiX, 
where X is any point in the two-dimensional space. Sim- 
ilarly, we define angular variables </>^ at A\. We also in- 
troduce 4>o(X) as the angle between axis OA\ and vec- 
tor OX. The solution in terms of these new variables is 
given by 



= 2(/) - 



1 



v 

£( 

i=l 



(A2) 



Equation (|A2|) satisfies Eq. §3§ because 9(x, y) = 
arctan ^-a-o ) is a solution of Laplace's equation for ar- 
bitrary (xo,yo). It also satisfies the strong anchoring 
boundary condition on the droplet due to a simple ge- 
ometrical fact (see Fig. IT2"]) : 

If C is a circle of radius a centered at 
point O, A is any point outside C and A' is 
the image of A under inversion, and B is any 
point on C, then ABO A = ABA' A - ZOAB. 



With this observation, we conclude that at the bound 
ary, <j>k+<j>'^ = 4>o+'K — 2ir(k — l)/p, Then 8 at the bound 
ary is given by 6 = 4>o + n/p — l/p[pn — n(p - 
i.e., the boundary conditions are satisfied. 

To calculate the free energy it is advantageous to in- 
troduce another field x, such that diX = e.ikdkO, [see 
Eq. JSJ], where eik is the antisymmetric tensor in two 
dimensions, e xx = e yy = 0, e xy = —£ yx = 1. It is easy to 



see that the free energy in terms of the Cauchy conjugate 
function \ has the same functional form as Eq. 



(A3) 
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j = p — a 



FIG. 12: Illustration of the geometrical fact that ABO A 
ABA' A - AOAB. 



The solution in terms of x is given by 



i p > > 

/(r) =21n|r| - - ^(ln |r - OAi\ + In |r - OA'A\) (A4) 



Upon integrating by parts in Eq. (|A3[) . we find 



(A5) 



where n is an outward unit normal on the circle and 
the contour C includes both the circle and a boundary 
at infinity. It is clear that for our solution the surface 
integral at infinity vanishes. 

Note that in polar coordinates around the center of the 
circle Mo), V0(r) = nf + ^^. Therefore V X (r) = 
j-§t^ — 0o ff, because by definition, V% is V# ro- 
tated by — 7r/2. Since on the boundary 6 = (f>o, we have 



/ x(V • X n)dl = / x*h- (A6) 

This equation can be further simplified with the help of 
the identity 



l7r In \u 2 + v 2 — 2uv cos < 
lo ln(max{|u|, \v\}) 

which leads to 



4tt, (A7) 



(A8) 



To evaluate the second term in cq. (|A5|) we 
use V 2 In |r - r | = 27r<5^ 2 ^(r - r ) and the fact that all 
defects contribute equally due to symmetry. Thus up to 
core energies that depend on microscopic details we find 



/ 



2tt 



2 \ 



xV 2 xd 2 x = - 4tt ln(R) + — In [R 

P \ RJ 

-£(ln|0X-0^| + 



(A9) 



InpA't-OAx 



We simplify this formula by taking advantage of com- 
plex numbers: let OAi = Rzi and OA- = ^-z,, 
where = exp[27ri(fc — l)/p]- Then we use the following 
chain of identities; valid for Z\ = \: 





- Zi) 


= z p - 1 =>, 












p-i 




-Zi) 




i=2 






l> 




p-1 


i=2 


- Zi) 


= E 1= ^ 

i=0 



(A10) 



to sum up the terms in the previous equation, which re- 
duces to 



- / xV 2 xd 2 x = - In 



p \p{R 2 P - a 2 P) 



pE c , (All) 



where we put the core energy E c back in and in- 
troduced the core radius c. Upon combining the 
Eqs. (|A11[) and (j A8|) we obtain the free energy ([5]), which 
upon minimization gives Eq. 



APPENDIX B: ORDERED DROPLETS IN 
DISORDERED FLUID 

In this appendix we briefly discuss the equilibrium con- 
figuration of defects inside a droplet of anisotropic two- 
dimensional liquid embedded in disordered host fluid. As 
before, we assume strong anchoring boundary conditions, 
which imply that the net charge of the defects in the in- 
terior must sum up to +1. Therefore we need p point 
defects with charges I /p. It is natural to assume that 
these can be arranged around the center of the droplet 
cither in regular p-gon or in regular (p — l)-gon with a 
defect in the center. We can then calculate preferred 
separation between the defects and the center of the cir- 
cle, R = Sa. This can be done as in Appendix [AJ and 
the results are given below. 
Regular p-gon, 
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ng ( (a?P - R 2 P)R?- 1 pc 

F — 111 



P 



which upon minimization gives 



>rl-> 



(Bl) 



5 = 



Regular (p — l)-gon (with central defect). 



ng(p-l) {{a 2 ?- 2 - R 2 P~ 2 )RP(p-r 
s In 



-,3p-2 



— In ' 



(B3) 



which leads to 



P 



3p-2 



(B4) 



Our calculation shows that the p-gon is energetically 
favorable compared to a (p — l)-gon for all values of p. 
Depending on anchoring energy, it is also possible to have 
surface and virtual defects instead of topological ones, 
but we do not pursue this possibility further here (see 
Ref. [3(| for an analysis when p = 1). 



APPENDIX C: SOLUTION FOR ^ FIELD 



the origin, where Laplace's equation can be solved eas- 
ily. Upon passing to complex coordinates £ and z, this 
transformation is 



(CI) 



(B2) where the real parameter w is given by 



r 2 + a\ — a 2 
' 2~r 



^{a 2 -a 2 ) 2 +r 2 [r 2 -2(a 2 + a 2 )} 



(C2) 



2r 



The outer radius of the ring is a\ and the inner radius a 
is 



i\ - w(r - a 2 ) 

r — a,2 — w 



(C3) 



The solution of Laplace's equation in this annular domain 
is then given by 



«2 H\C\/ai)~aiH\C\/a) 
ln(ai/a) 



(C4) 



which can be easily expressed in terms of z — (x, y) by 
inverting Eq. ([CI]) . Now, Vj, = ^J(Vtp) 2 d 2 x can be 
calculated as follows: 



In this appendix we calculate the droplet-droplet po- 
tential contribution We assume that droplets of 
radii a\ and 02 are located at (0, 0) and (r, 0) respec- 
tively, and that ip satisfies Laplace's equation everywhere 
in space, rp is forced to equal to — a± and — ct2 on the 
surface of the droplets. This problem can be solved by 
taking advantage of the fact that conformal mappings 
leave Laplace's equation invariant. Hence, we search for 
a conformal mapping that transforms two circles (repre- 
senting the droplet boundaries) into an annular ring at 



- 9 -f Vi>.Ml-if 

* J droplet one * J droplet two 



' droplet one 

■ng{a\ - a 2 ) 2 
ln(ai/a) 



Vij) ■ Ml 
V 2 ipd 2 x 



droplet two 



(C5) 



In the limit r>Oi, a 2 this equation reduces to Eq. (fl 
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